AD-A079  5*3  AEROSPACE  CORP  EL  SE9UND0  CA  F/6  17/7 

COMBINED  PSEUDO  RAN9E  AND  DOPPLER  P0SITI0NIN8  FOR  THE  STATIONAR— ETCIU) 
MAR  79  P  S  J0R8ENSEN  F0A701-78-C-0079 


|«NCLASSIFIEO  TR-0079<  4975-01) -3  SAMSO-TR-79-79  NL 


END 


DATE 

FilMEO 

2-80 


ooc 


MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BURtAU  Of  STANDAROft-1963-.1 


®A079543 


Combined  Pseudo  Range  and 
Doppler  Positioning  for  the 
Stationary  Navstar  User 


Prepared  by 

P.  S.  JORGENSEN 
Programs  Group 
The  Aerospace  Corporation 
El  Segundo,  Calif.  90245 


20  March  1979 


Final  Report 


\n\ 


.:,r*  i 


- 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


Prepared  for 


SPACE  AND  MISSILE  SYSTEMS  ORGANIZATION 
AIR  FORCE  SYSTEMS  COMMAND 
Los  Angeles  Air  Force  Station 
P.0,  Box  92960,  Worldway  Postal  Center 
Los  Angeles.  Calif.  90009 


80 


1 


o  ^ 


This  final  report  was  submitted  by  The  Aerospace  Corporation, 

El  Segimdo,  CA  90245,  under  Cortract  F04701-78-C-0079  with  the  Space 
and  Missile  Systems  Organization,  Deputy  for  Space  Communications 
Systems,  P.O.  Box  92960,  Worldway  Postal  Center,  Los  Angeles,  CA  90009. 
It  was  reviewed  and  approved  for  The  Aerospace  Corporation  by  E.  M. 
Lassiter.  Programs  Group. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  (PAS)  and  is 
releasable  to  the  National  Techi'-ical  Information  Service  (NTIS).  At  NTIS 
it  will  be  available  to  the  general  public,  including  foreign  nations.  Pub¬ 
lication  of  this  report  does  not  constitute  Air  Force  approval  of  the  report's 

findings  or  conclusions.  It  is  published  only  for  the  exchange  and  stimulation 
of  ideas. 

^  y;  // 


,,^SEPH  A.  STRAOA,  LCDR,  USN 
/  Director  of  Navioatlon  Eouipment  &  Avionics 


FOR  THE  COMMANDER 


J^PH  C.  CARLING,  Lt  Col,  Uw 

Assistant  Deputy  for  Space  Navigation  Systems 


U\C  LASSIFIED 


security  classification  of  this  page  ClWlM  Data  emmrmd) 


REPORT  DOCUMENTATION  PAGE 


REPORT  NUM  0  E  R/  I 

SAMS(^-|tV^  9-7  9 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


/ 


2.  GOVT  ACCESSION  NO 


3.  RECIPIENT'S  CATALOG  NUMBER 


Title  fand  Subtitle) 

\  Combined  Pseudo  Range  and  Doppler 

Positioning  for  the  Stationary  Navstar  Use 


S._TVPE  OF  REPORT  A  PERIOD  COVERED 


Final 


PERFORMIMC  ORG-  REPORT  NIIURF^B 

TR-O079(4475-^1  )-3 


7.  authoR(s; 

Paul  S,/  Jorgensen 


ibj  p'-'-'s./ 


lllfL: 


F047Ol-78-C-J0^79 


PERFORMING  organization  NAME  AND  ADDRESS 

The  Aerospace  Corporation 
El  Segundo,  Calif.  90245 


0.  program  element,  PROJECT,  TASK 
AREA  a  WORK  UNIT  NUMBERS 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Space  and  Missile  Systems  Organization 

Air  Force  Systems  Command 

Los  Angeles  Air  Force  Station 

P.  O.  Box  Q2960 

Los  Angeles,  Calif.  90009 


2-  REPORT  PATE 

[ar|MII7  9 


IS.  security  class.  (oI  Ihia  rtporl) 
Unclassified 


l«  MONITORING  AGENCY  NAME  A  AOORESSri/ <JI//aran(  Irom  CentrolllnS  Ottic*) 


ISa.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


IS.  DISTRIBUTION  ST  ATEMEN  T  fof  fh/a  RaporO 


Approved  for  public  release;  distribution  unlimited. 


t7.  distribution  statement  (of  the  ebetfsct  entered  in  Block  20,  if  different  from  Report) 


ta.  supplementary  notes 


19.  KEY  WORDS  ('Continue  on  reveree  aide  if  neceeaery  and  Identify  by  bloclr  numboO 

Biphase  Modulation  Manpack  Receiver 

Doppler  Positioning  Navstar 

Global  Positioning  System  Pseudorandom  Noise  Modulation 

Integrated  Doppler  PRN  Code  Tracking 

L-band  Carrier  Phase  Lock  Loop  Satellite  Geodesy 


20  ABSTRACT  (Continue  on  ravarae  a/<fa  H  naeaaaaiy  ond  Idonllly  by  block  number) 

'An  early  success  in  testing  Navstar  user  equipment  was  achieved  with  the  Man- 
pack  receiver.  Using  only  two  satellites,  Navstars  1  and  2,  the  receiver  accu¬ 
rately  performed  two-dimensional  positioning  while  stationary.  It  was,  in 
effect,  combining  pseudo  range  and  integrated  Doppler  data.  This  concept  of 
combining  both  types  of  data  can  be  extended  to  three  dimensions  and  to  the 
one  rational  Phase  III  Navstar.  The  stationary  user  can  accomplish  three-  - 


•v/ 


DO  U73 

rrACSiMiLCi 


UNCLASSIFIED  OALJ'OO 


security  classification  of  This  rage  firh»n  note  Entered)  ) 


T 


UNCLASSIFIED 


SECURITY  CLASSIFICATIOH  OF  THIS  PAGErniFn  Dmtm  Entmna) 


Abstract  (Continued) 


dimensional  navigation  with  only  three  satellites.  This  capability  is  of  inter¬ 
est  to  the  Army  since  it  would  permit  the  stationary  use  of  Navstar  with  much 
greater  elevation  mask  angles.  With  the  baseline  Z4-satellite  constellation, 
three-dimensional  navigation  can  be  accomplished  with  a  30-deg  mask  angle,  but 
it  will  take  a  few  minutes  longer  because  of  the  extra  time  required  for  making 
integrated  Doppler  measurements.  Using  this  technique,  three-dimensional 
position  will  be  obtained  to  within  16  meters  about  90  percent  of  the  time.  This 
operational  capability  will  enhance  the  Army's  ability  to  use  Navstar  in  the 
presence  of  jamming,  foliage,  and  terrain. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  RAOErWAwi  Dmim  Enlmnd) 


:»aiaui 


CONTENTS 


I,  INTRODUCTION . 

II.  DISCUSSION  AND  RESULTS  . . 

in.  ACCURACY  . 

IV.  COMPUTATIONAL  TECHNIQUES  . 

V.  EQUATIONS  USING  BOTH  PSEUDO  RANGE 

AND  INTEGRATED  DOPPLER  . 


APPENDICES; 

A.  SIGNAL  STRUCTURE . 

B.  NAVIGATION  EQUATIONS  . 

C.  BASELINE  ORBIT  CONFIGURATION  .. 

D.  DOPPLER  NAVIGATION  WITH  NAVSTAR 


1 


I.  INTRODUCTION 


The  Navs tar /Global  Positioning  System  (GPS)  is  a  satellite -based 
navigation  system  that  will  provide  extremely  accurate  three-dimensional 
position  and  velocity  information  to  properly  equipped  users  anywhere  on  or 
near  the  earth.  It  is  a  Joint  Service  Program  managed  by  the  USAF  with 
deputies  from  the  Navy,  Army,  Marines,  Defense  Mapping  Agency,  and 
Coast  Guard,  and  with  technical  support  provided  by  The  Aerospace 
Corporation.  The  baseline  program  is  divided  into  three  phases: 

I  -  Concept  Validation  Phase  (1974-1979) 

II  -  Full-Scale  Engineering  Development  Phase  (1979-1983) 

LII  -  Production  Phase  (1983-  ) 

The  current  baseline  orbital  configuration  for  the  fully  operational 
Phase  III  GPS  employs  24  satellites  in  55-deg  inclined,  circular,  12-hr 
orbits  to  transmit  navigation  signals.  Continuous  three-dimensional  global 
coverage  will  be  provided  by  placing  eight  satellites,  equally  spaced,  in  each 
of  three  orbit  planes  120  deg  apart  in  inertial  space  (Fig.  1). 

The  major  elements  comprising  the  navigation  payload  on  the  satellites 
are  the  pseudorandom  noise  signal  assembly  (PRNSA),  atomic  frequency 
standard,  processor,  and  Li-band  antenna.  The  PRNSA  includes  the  baseband 
generator,  which  produces  the  basic  P  (precise)  and  C/A  (coarse /acquisition) 
ranging  codes  and  encodes  navigation  data  from  the  processor  onto  the  pseudo¬ 
random  noise  (PRN)  ranging  signal;  the  amplifier /modulator  units,  which 
supply  the  Li  (1575  MHz)  and  Lz  (1227  MHz)  carrier  frequencies  modulated 
by  the  PRN  ranging  signals;  and  the  high  power  amplifiers,  which  amplify 
the  carrier  signals  for  transmission. 

The  Control  Segment  consists  of  a  Master  Control  Station  (MCS),  a 
navigation  message  Upload  Station  (ULS),  and  widely  separated  Monitor 


SPACECRAFT  CONFIGURATION 

•  950  lb  (typicil) 

•  300  W  (typicili 

•  NAV  SIGNALS 

•  1227  MHz 

•  1575  MHz 

•  ONE-WAY  TRACKING 

•  lO'’^  CLOCK  STABILITY 


Fig.  1.  Orbital  Configuration 


4 


Stations  (MSs).  The  MSs  passively  track  all  satellites  in  view  and  accumu¬ 
late  ranging  data,  which  is  processed  at  the  MCS  to  calculate  the  satellite 
ephemerides,  clock  drifts,  and  propagation  delay.  At  least  once  a  day  this 
information  is  transmitted  by  the  ULS  to  the  satellites  for  subsequent  down¬ 
link  transmission  of  the  navigation  data  encoded  on  the  carrier  signals. 

The  conventional  user  measures  pseudo  range  and  pseudo  range  rate 
using  the  navigation  signal  from  each  of  four  satellites.  "Pseudo  range" 
is  the  true  distance  from  the  satellite  to  the  user  plus  an  offset  due  to  the 
user's  clock  bias.  Similarly,  "pseudo  range  rate"  is  the  true  slant  range 
rate  plus  an  offset  due  to  the  frequency  of  the  user's  clock.  Each  signal 
carries  ephemeris  data  and  system  timing  information  for  that  satellite. 
This  allov/s  the  user  receiver /proce  ssor  to  convert  the  pseudo  range  and 
pseudo  range  rate  to  user  three-dimensional  position  and  velocity. 

Obviously,  a  stationary  user  has  no  need  to  determine  his  velocity,  so 
only  pseudo  range  data  is  needed.  However,  pseudo  range  rate  (Doppler 
frequency  offset)  provides  additional  position  information.  The  stationary 
user  has  the  special  opportunity  to  simultaneously  process  both  types  of  data 
to  determine  his  position,  which  allows  position  determination  in  three 
dimensions  using  only  three  satellites  instead  of  four.  To  obtain  sufficiently 
accurate  integrated  Doppler  measurements,  the  stationary  user  must  op¬ 
erate  his  receiver  a  few  minutes  longer. 
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II.  DISCUSSION  AND  RESULTS 


By  combining  pseudo  range  and  integrated  Doppler  data,  the  stationary 
user  has  the  advantage  of  two  types  of  navigation  data  even  though  only  one 
satellite  is  tracked.  From  the  pseudo  range  data  the  user  obtains  position 
information  just  as  the  conventional  Navstar  user  does,  and  from  the 
integrated  Doppler  measurements  he  obtains  position  information  via  the 
technique  geodetic  users  employ  with  Transit  satellites.  Each  Navstar 
satellite  serves  as  two,  allowing  for  both  the  usual  Navstar  pseudo  ranging 
and  Transit-like  Doppler  positioning.  This  duality  of  operation  is  particu¬ 
larly  effective  because  the  "two  satellites"  appear  to  the  user  to  be  90  deg 
apart.  The  pseudo  range  measurement  provides  position  information  along 
the  line  of  sight  (LOS)  to  the  satellite,  and  the  integrated  Doppler  measure¬ 
ment  provides  information  perpendicular  to  this  LOS  (it  provides  position 
information  along  the  direction  of  the  component  of  relative  velocity  that  is 
perpendicular  to  the  LOS). 

The  conventional  Navstar  user  requires  pseudo  range  measurements 
from  four  satellites  because  he  must  compensate  for  the  user  clock  phase 
bias,  which  constitutes  a  fourth  unknown  in  the  navigation  computations  (in 
addition  to  the  three  oosition  components).  The  stationary  user  using  both 
types  of  data  is  confronted  with  a  fifth  unknown  in  his  navigation  solution: 
the  user  clock  frequency  bias.  For  this  reason  three  satellites  must  be 
used  to  obtain  at  lease  five  independent  measurements. 

The  stationary  user  could  perform  pseudo  ranging  to  three  satellites 
and  integrated  Doppler  measurements  to  only  two.  An  alternate  approach 
is  to  make  both  types  of  measurement  to  all  three  satellites  and  from  these 
six  independent  measurements  perform  a  least  squares  solution  of  the  five 
unknown  quantities.  This  latter  approach  is  the  method  used  in  generating 
the  quantitative  results  given  in  this  report. 
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One  interesting  result  of  these  computations  is  that  whenever  there 
are  only  three  satellites  visible  to  the  stationary  user  (because  of  either  a 
very  high  elevation  mask  angle  or  a  small  number  of  satellites  in  the 
constellation),  the  user  always  obtains  a  good  solution.  This  differs  from 
the  conventional  approach  in  which  the  user  may  not  get  a  good  navigation 
solution  when  only  four  satellites  are  available  because  the  four  satellites 
as  a  group  have  a  poor  geometry  relative  to  the  user. 

The  power  of  using  both  types  of  data  arises  from  two  factors.  First, 
the  navigation  solutions  (as  simulated  for  this  study)  are  obtained  from  an 
overdetermined  set  of  equations  (six  equations  with  five  unknowns).  Second, 
the  pseudo  range  and  integrated  Doppler  data  are  inherently  complementary. 
This  can  be  illustrated  by  considering  a  seemingly  poor  geometric  arrange¬ 
ment  of  only  three  satellites,  all  at  rather  high  elevation  angles  relative 
to  the  user.  Let  us  consider  the  case  in  which  one  satellite  crosses  the 
orbit  plane  of  the  other  satellites  in  the  constellation  and  these  other  two 
satellites  are  on  opposite  sides  of  the  first  (crossing)  satellite:  These  three 
satellites  are  lined  up  in  a  row,  so  their  geometry  does  indeed  seem  quite 
poor  to  a  user  directly  beneath  them.  However,  the  velocity  of  the  middle 
satellite  is  almost  perpendicular  to  that  of  the  other  two  satellites  and  is 
also  perpendicular  to  the  plane  defined  by  the  instantaneous  position  of  all 
three.  This  direction  of  the  velocity  of  the  middle  satellite  is  precisely 
the  additional  information  needed  by  the  stationary  user  to  obtain  a  position 
component  perpendicular  to  the  plane  common  to  himself  and  the  three 
satellites . 

The  stationary  user  obtains  a  good  solution  as  long  as  three  satellites 
are  available.  This  is  illustrated  by  three  cases:  In  the  first  case  three 
satellites  within  the  baseline  24-satellite  constellation  are  visible  above  an 
elevation  mask  angle  of  30  deg  about  99.8  percent  of  the  time,  and  the  user 
does  indeed  obtain  a  good  navigation  solution  99.  8  percent  of  the  time.  The 
cumulative  probability  of  position  dilution  of  precision  (PDOP)  is  shown  in 
Fig.  2,  (The  meaning  of  PDOP  is  discussed  in  detail  later;  for  the  present 
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consider  it  to  be  an  error  ratio,  the  ratio  of  the  user's  three-dimensional 
position  error  to  his  error  in  ranging  to  the  satellite.  )  The  PDOP  generally 
lies  between  2.5  (10th  percentile)  and  3,  1  (90th  percentile).  Since  the 
ranging  accuracy  of  the  overall  Navstar  system  is  expected  to  be  about  5.3  m 
the  stationary  user  will  obtain  his  three-dimensional  position  to  within  16  m 
90  percent  of  the  time. 

In  the  second  case,  with  only  12  satellites  and  a  5-deg  elevation  mask 
angle,  three  satellites  are  visible  about  97  percent  of  the  time,  and  once 
again  the  stationary  user  obtains  a  good  navigation  solution  for  all  cases 
when  at  least  three  satellites  are  visible  (Fig.  3  shows  the  PDOP  statistics 
for  this  case).  Again,  the  user  obtains  his  position  to  within  16  m  90  per¬ 
cent  of  the  time. 

In  the  third  case,  with  15  satellites  and  a  10-deg  elevation  mask  angle, 
three  satellites  are  visible  99  percent  of  the  time,  and  the  user  obtains  his 
position  to  within  14  m  90  percent  of  the  time  (the  PDOP  statistics  for  this 
case  are  shown  in  Fig.  4). 


Fig  2.  PDOP  for  24-Satellite  Constellation 
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III.  ACCURACY 


Pseudo  range  measurements  are  made  by  tracking  the  PRN  biphase 
modulation  of  the  L-band  carrier  (see  the  discussion  of  signal  structure  in 
Appendix  A).  The  chipping  rate  of  this  modulation  is  10.23  MHz,  corre¬ 
sponding  to  a  wavelength  of  about  30  m.  The  integrated  Doppler  is  measured 
with  a  phase  lock  loop  that  tracks  the  Lj  carrier  whose  frequency  is  1575.42 
MHz,  which  corresponds  to  a  wavelength  of  about  19  cm.  The  granularity 
of  these  two  tracking  operations  thus  differs  by  about  two  orders  of  magni¬ 
tude,  and  the  accuracy  of  the  integrated  Doppler  can  therefore  be  expected 
to  be  much  greater  than  that  obtained  with  pseudo  ranging.  This  high 
degree  of  integrated  Doppler  measurement  accuracy  is  crucial  to  the  sta¬ 
tionary  Navstar  user  who  wishes  to  obtain  an  accurate  position  fix  in  a 
reasonable  amount  of  time  via  both  types  of  measurements.  (Incidentally, 
this  highly  accurate  carrier  tracking  is  also  important  to  current  geodetic 
users,  who  now  obtain  accurate  position  fixes  using  Transit  satellites.  ) 

The  analysis  of  satellite  constellation  performance  with  a  conventional 
Navstar  user  (who  uses  pseudo  range  measurements  to  four  satellites)  uses  a 
technique  referred  to  as  geometric  dilution  of  precision  (GDOP),  which  is  dis¬ 
cussed  in  Appendix  B.  The  quantity  generally  of  interest  is  PDOP,  which  is 
a  dimensionless  error  ratio  (the  ratio  of  three-dimensional  user  position 
error  to  pseudo  ranging  error).  The  baseline  (24-satellite)  constellation  is 
defined  in  Appendix  C,  which  includes  a  graph  showing  the  statistical  cumu¬ 
lative  probability  of  PDOP  for  this  constellation. 

For  the  case  of  a  stationary  user  who  uses  both  pseudo  range  and 
integrated  Doppler  measurements,  it  would  be  convenient  to  define  a 
similar  nondimensional  error  ratio  to  describe  performance.  This  is 
complicated  by  the  fact  that  we  are  now  dealing  with  two  types  of  measure¬ 
ments  that  may  not  be  of  comparable  accuracy.  In  this  report  it  was  assumed 
that  the  integrated  Doppler  measurements  are  50  times  more  accurate  than 
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the  pseudo  range  measurements.  (For  example,  if  the  accuracy  ( lo)  in 
pseudo  ranging  were  about  5.3  m.  the  integrated  Doppler  measurement 
would  be  accurate  to  within  about  11  cm.  This  level  of  Doppler  performance 
is  consistent  with  current  Navstar  test  results.  )  The  PDOP  is  thus  defined 
as  the  ratio  of  three-dimensional  navigation  error  to  pseudo  ranging  error, 
with  the  integrated  Doppler  error  considered  to  be  one -fiftieth  that  of  the 
pseudo  range  measurement  errors. 

Another  consideration  is  the  length  of  time  over  which  the  integrated 
Doppler  measurement  should  be  taken  (about  2  min  seems  appropriate; 
during  this  time  the  satellites  will  have  rotated  in  their  orbits  about  1  deg). 

It  has  been  shown  that  using  this  2-min  assumption  and  the  50:1  accuracy 
ratio  assumption  leads  to  a  solution  in  which  the  errors  in  pseudo  range 
and  integrated  Doppler  measurements  each  contribute  about  equally  to  the 
stationary  user  three-dimensional  navigation  error. 

The  quantitative  results  given  in  Sec.  II  are  based  on  computations  in 
which  it  is  assumed  that  three  satellites  are  simultaneously  tracked  for 
2  min.  (Stationary  users  generally  have  sequential  receivers,  so  three 
integrated  Doppler  measurements  require  6  min.)  Two  points  about  the 
"extra"  time  the  stationary  user  needs  to  perform  integrated  Doppler 
measurements  should  be  noted.  First,  in  many  cases  only  three  pseudo 
range  and  two  Doppler  measurements  are  strictly  necessary,  which  would 
increase  the  time  by  only  4  min.  Second,  during  this  "extra"  time  the 
stationary  user's  receiver  is  oerforming  other  necessary  functions,  e.g., 
reading  the  navigation  data  message. 

Note  that  position  could  be  obtained  using  only  integrated  Doppler 
measurements  from  four  satellites;  this  parallels  the  conventional  case 
of  four  pseudo  range  measurements.  Although  it  is  outside  the  scope  of 
this  report,  this  approach  has  been  investigated  and  is  discussed  in 
Appendix  D. 


IV.  COMPUTATIONAL  TECHNIQUES 


A  computer  program  has  been  developed  that  divides  the  earth  into 
areas  bounded  by  increments  of  latitude  and  longitude  and  computes  the 
GDOP  parameters  for  a  user  at  the  center  of  these  areas.  The  program 
determines  which  satellites  are  visible  to  the  user  and  searches  through  all 
combinations  to  find  the  set  of  three  satellites  that  yields  the  best  GDOP 
parameters.  In  general,  the  criterion  most  often  used  is  the  three-dimen¬ 
sional  PDOP.  To  take  into  account  satellite  motion  and  earth  rotation,  the 
computations  are  repeated  for  all  latitudes  and  longitudes  at  successive 
time  increments;  thus,  the  statistics  of  the  GDOP  parameters  are  obtained 
for  any  time  of  day  and  any  place  on  earth.  In  general  the  statistics  are 
computed  on  a  worldwide  basis,  giving  equal  weight  to  all  locations  on 
earth.  These  statistics  are  used  to  obtain  cumulative  percentiles. 

For  each  of  the  statistical  results  given  in  this  report,  a  total  of 
3,  Z40  computations  is  performed.  For  each  of  these,  the  best  three 
satellites  are  obtained,  and  the  PDOP  is  computed.  The  Northern  Hemis¬ 
phere  is  divided  into  18  latitudes,  separated  by  5  deg,  and  12  longitudes, 
also  separated  by  5  deg.  The  computations,  which  are  performed  for  the 
center  of  each  of  these  'boxes,  '  are  repeated  15  times  over  a  time  span 
during  which  the  satellites  rotate  through  an  angle  equal  to  the  spacing 
between  satellites  in  a  common  orbit  plane.  This  effectively  covers  all 
locations  on  the  earth's  surface  and  all  times  of  day. 
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V.  EQUATIONS  USING  BOTH  PSEUDO  RANGE 
AND  INTEGRATED  DOPPLER 


The  conventional  Navstar  technique  of  navigating  by  means  of  pseudo 
ranging  to  four  satellites  is  discussed  in  Appendices  B  and  C.  Equations 
that  apply  to  the  case  in  which  both  pseudo  range  and  integrated  Doppler 
measurements  are  employed  are  presented  here. 

Making  pseudo  range  measurements  at  the  beginning  of  a  time  interval 
and  performing  integrated  Doppler  measurements  during  this  interval  is 
mathematically  equivalent  to  making  two  sets  of  pseudo  range  measurements, 
one  set  at  the  beginning  and  the  other  at  the  end  of  the  interval.  The  inte¬ 
grated  Doppler  measurements  are  the  same  as  the  difference  between  the 
two  sets  of  pseudo  range  measurements  to  the  same  satellite.  (In  theory, 
one  could  imagine  an  approach  using  only  pseudo  range  measurements, 
enabling  the  user  to  employ  fewer  than  four  satellites.  The  difficulty  with 
this  approach  lies  in  the  random  errors  from  the  pseudo  range  measurements; 
for  this  approach  to  work  with  reasonable  accuracy,  the  time  between  the 
two  sets  of  measurements  would  have  to  be  about  an  hour.  )  Using  integrated 
Doppler  measurements,  which  take  advantage  of  the  very  high  accuracy 
with  which  the  L-band  carrier  can  be  tracked,  the  user  can  obtain  an  accu¬ 
rate  position  determination  in  just  a  few  minutes.  However,  from  the  point 
of  view  of  developing  the  equations  associated  with  combining  the  two  types 
of  measurements,  this  process  can  be  considered  the  mathematical  equiva¬ 
lent  of  two  sets  of  pseudo  range  measurements. 

Linearized  solution  equations  for  obtaining  user  position  from  pseudo 
range  measurements  only  have  been  developed  and  conveniently  expressed 
in  matrix  notation  (Appendix  B).  The  case  of  two  time  points  and  three 
satellites  involves  six  equations  in  five  unknowns  that  relate  user  position 
and  clock  offsets  to  pseudo  range  measurements  are  given  as 
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where  is  the  direction  cosine  of  the  angle  between  the  range  to  the  ith 
satellite  and  the  jth  coordinate  at  the  kth  time;  Ax,  Ay,  and  Az  are  the 
differences  between  the  actual  and  nominal  (a  priori  best  estimate)  values  of 
user  position;  ATj^  is  the  difference  between  the  actual  and  nominal  values 
of  user  clock  bias  at  the  kth  time  (AT^  and  AT^  will  differ  because  of  the 
offset  of  the  user's  clock  frequency  from  nominal);  and  is  the  difference 

between  the  actual  and  nominal  pseudo  range  measurement  to  the  ith  satellite 
at  the  kth  time . 


The  last  three  pseudo  range  equations  at  the  second  time  point  can  be 
replaced  with  integrated  Doppler  equations  by  taking  the  difference  between 
each  of  these  equations  and  the  corresponding  equation  at  the  first  time 
point.  Let 
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where  is  the  change  in  the  direction  cosine  of  the  angle  between  the 
range  to  the  ith  satellite  and  the  jth  coordinate  over  the  time  interval. 
At  is  the  change  in  the  difference  between  the  actual  and  nominal  values 
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of  user  clock  bias  over  the  time  interval  (the  average  offset  of  the  user's 
clock  frequency  from  nominal  multiplied  by  the  time  interval),  and  is 
the  change  in  the  difference  between  the  actual  and  nominal  pseudo  range 
measurement  to  the  ith  satellite  over  the  time  interval  (the  difference 
between  the  actual  and  nominal  integrated  Doppler  measurements). 

With  these  substitutions,  the  six  equations  in  five  unknowns  that 
relate  user  position,  clock  bias,  and  average  clock  frequency  offset  to 
pseudo  range  and  integrated  Doppler  measurements  are  given  by 
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Using  matrix  notation,  the  above  equations  can  be  expressed  com¬ 
pactly.  Let 

r  =  the  six -element  pseudo  range  and  integrated  Doppler 
measurement  difference  vector 

X  =  the  user  position  and  clock  correction  vector 

A  =  the  6x5  solution  matrix 


19 


^111 

'^121 

‘^131 

1 

0 

“211 

"221 

“231 

1 

0 

A  s 

®311 

*^321 

°'331 

1 

0 

^12 

®13 

0 

1 

®21 

^22 

®23 

0 

1 

®31 

^32 

S3  3 

0 

1 

X  = 

|~Ax 

Ay 

Az 

ATf 

At 

r  = 

[}«.i 

AR21 

ADj 

AD 

->  AD~| 
2 

Therefore 

Ax 

=  r 

which  compactly  expresses  the  relationship  between  pseudo  range  and 
integrated  Doppler  measurements  and  user  position,  clock  bias,  and  clock 
frequency  offset.  Since  this  relationship  is  linear,  it  can  be  used  to  express 
the  relationship  between  the  errors  in  the  measurements  and  the  user  quan¬ 
tities  and  can  be  stated  as 


£ 

r 


where  represents  the  measurement  errors  and  e  the  corresponding  errors 
in  user  position,  clock  bias,  and  clock  frequency  offset. 

Let  us  now  consider  the  covariance  matrices  of  the  expected  measure¬ 
ment  errors  and  the  user  quantities.  The  first  covariance  measurement  is 
a  6  X  6  array  composed  of  the  expected  values  of  the  squares  and  products 
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of  the  measurement  errors.  The  diagonal  terms  in  the  matrix  (the  squares 
of  the  expected  errors)  are  the  variances  (i.e,,  the  squares  of  the  expected 
one  sigma  measurement  error  values).  The  off-diagonal  terms  represent 
the  covariance  between  the  measurements  and  reflect  the  correlations  to 
be  expected.  Similarly,  the  covariance  matrix  for  the  user  quantities  is 
composed  of  the  expected  values  of  the  squares  and  products  of  the  errors 
in  the  user  quantities.  The  diagonal  terms  are  the  variance  or  the  squares 
of  the  errors  in  user  position  and  clock  parameters  (la),  while  the  off- 
diagonal  terms  reflect  the  correlations  in  these  errors.  These  covariance 
matrices  are 


COV(r)  =  E 
COV(x)  =  E  )  e  e 

I  ""  ) 

where  the  symbol  E  |  |  designates  the  expected  value  of  the  quantity  inside 

the  braces. 


When  a  weighted  least  squares  solution  is  performed  for  the  five  user 
quantities  from  the  six  measurements,  the  relationship  between  the  measure¬ 
ment  covariance  matrix  and  the  covariance  matrix  for  the  user  quantities 
can  be  shown  to  be 


COV(x)  =  [A^COV(r)'^Aj'‘ 

Note  that  the  relationship  between  the  measurement  errors  and  the  user's 
errors  is  a  function  only  of  the  solution  matrix  A,  which  in  turn  is  a  function 
only  of  the  direction  cosines  of  the  LOSs  from  the  user  to  the  satellites 
along  whichever  coordinate  system  is  used.  In  other  words,  the  error 
relationships  are  functions  only  of  satellite  geometry.  An  important  consid¬ 
eration  in  the  proper  use  of  Navstar  is  the  fact  that  the  three  satellites 
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being  used  must  possess  'good''  geometric  properties.  (A  'good  '  situation 
is  one  in  which,  because  of  satellite  geometry,  a  given  level  of  error  in  the 
measurements  results  in  small  user  errors.  ) 


Certain  assumptions  regarding  pseudo  range  measurement  errors 
lead  to  a  method  of  quantitatively  determining  whether  a  particular  three - 
satellite  geometry  is  good  or  bad.  Let  each  individual  pseudo  range  mea¬ 
surement  have  an  error  (la)  of  unity,  where  the  expected  mean  is  zero  and 
the  correlation  of  errors  between  satellites  is  also  zero.  If  the  integrated 
Doppler  measurements  are  one-fiftieth  those  of  the  pseudo  range  measure¬ 
ments  (as  discussed  previously),  the  reciprocal  of  the  measurement  covari¬ 
ance  matrix  is  given  by 
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Let  Vx,  Vy,  and  Vz  be  the  variances  of  the  three  components  of  user 
position.  The  three-dimensional  PDOP  is  then  given  by 

PDOP  =  \J  Vx  +  Vy  +  Vz 

The  "best"  set  of  three  satellites  is  that  set,  out  of  all  combinations 
visible  to  the  user,  that  has  the  minimum  PDOP.  The  cumulative  per¬ 
centiles  of  PDOP  for  various  combinations  of  total  number  of  satellites  and 
assumed  elevation  mask  angle  are  shown  in  Figs.  2  through  4. 
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APPENDIX  A 


SIGNAL  STRUCTURE 


Each  satellite  transmits  a  navigation  signal  on  two  L-band  frequencies, 
one  at  1575.42  MHz  (Lj)  and  the  other  at  1227.6  MHz  (U^).  These  two 
carrier  frequencies  are  biphase  modulated  by  pseudo  random  sequences 
providing  a  spread  spectrum  modulation.  The  carrier  is  actually  modu¬ 
lated  by  two  such  sequences  in  phase  quadrature  so  that,  strictly  speaking, 
this  carrier  is  actually  quadraphase  modulated.  One  pseudo  random  se¬ 
quence  is  a  precision  (P)  signal  at  a  random  pulse  repetition  rate  of  10.23 
MHz  and  is  an  extremely  long  code  so  that  for  all  practical  purposes  it  is 
a  truly  random  sequence.  The  second  pseudo  random  sequence  is  a  coarse/ 
acquisition  (C/A)  signal,  which  is  a  short  sequence  used  either  for  initial 
acquisition  of  the  P  signal  or  as  a  less  accurate  navigation  signal  for  low- 
cost  users.  The  carrier  frequency  is  biphase  modulated  only  by  the  P 
signal  or,  as  a  ground -controlled  option,  only  by  the  C/A  signal. 

For  the  purpose  of  this  discussion  we  will  confine  our  attention  only 
to  the  P  signal  on  the  carrier;  this  is  the  primary  navigation  signal. 

Note  that  the  carrier  frequency  ( 1  57 5.  42  MHz )  is  an  exact  multiple  (154)  of 
the  pseudo  random  sequence  pulse  rate  (10.23  MHz).  The  wavelength  of 
the  carrier  is  only  19  cm,  whereas  the  "chipping"  rate  of  the  pseudo  random 
sequence  is  about  30  m. 

The  pseudo  random  sequence  is  generated  by  a  feedback  shift  register, 
the  output  of  which  modulates  the  carrier  as  illustrated  in  Fig.  A-1.  The 
P  signal  pseudo  random  sequence  generator  is  functionally  illustrated  in 
Fig.  A -2.  By  combining  four  12-stage  feedback  shift  registers,  the 
equivalent  of  a  48-stage  shift  register  is  obtained.  Using  a  pseudo  random 
sequence  to  biphase  modulate  the  carrier  results  in  the  transmitted  spec¬ 
trum  being  spread  as  illustrated  in  Fig.  A-3.  This  spread  permits  an 
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interference  signal  from  being  rejected  in  the  receiver  in  the  following 
manner:  The  receiver's  pseudo  random  sequence  generator  modulates  the 
incoming  signal  in  the  same  manner  as  the  generator  on  the  satellite  (the 
two  generators  are  the  same).  The  original  carrier  frequency  is  thus 
reconstructed  and  is  collapsed  to  a  very  narrow  band.  The  interference 
signal,  however,  is  spread  out  over  a  wide  spectrum,  and  this  signal  can 
therefore  be  filtered  out  so  that  only  a  small  residue  remains  near  the 
now -reconstructed  carrier  frequency. 


m 


f 


m 


Fig.  A-1.  Pseudo  Random  (Noise)  Code 
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CHARACTERISTICS 
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CORRELATION 

PROPERTIES 


CODE  PROPERTIES 
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Fig.  A-2.  Code  Generator 
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Fig.  A-3.  Pseudo  Random  Noise -Spread  Spectrum 


25 


Functionally,  Navstar  receivers  incorporate  two  tracking  loops  that 
must  operate  simultaneously  to  properly  track  the  Navstar  navigational 
signal.  The  first  is  the  code  tracking  loop,  which  tracks  the  pseudo  random 
sequence  by  matching  the  locally  generated  sequence  with  the  sequence  on 
the  received  signal.  Simultaneously,  a  phase  lock  loop  is  tracking  the 
carrier  frequency.  The  actual  process  is  much  more  complicated  because 
of  the  necessary  intermediate  frequency  (IF)  downconversion  steps.  For 
simplicity,  however,  these  downconversion  steps  are  omitted  in  the  func¬ 
tional  diagram  shown  in  Fig.  A-4.  Figure  A -5  expands  on  the  function  of 
such  a  receiver  by  illustrating  the  role  of  a  feedback  shift  register  and 
envelope  detectors  in  the  code  lock  loop  that  tracks  the  incoming  pseudo 
random  sequence.  There  are  three  outputs  of  the  feedback  shift  register; 
an  on-time  sequence  P^,  an  early  sequence  P^,  and  a  late  sequence  P^. 

The  early  and  late  codes  modulate  the  carrier  frequency  C,  which  is  syn- 
thesizeu  in  the  phase  lock  loop.  These  signals  are  then  mixed  with  the 
incoming  signal,  thereby  generating  voltages  proportional  to  the  extent 
that  the  sequences  match  the  incoming  sequence  from  the  satellite.  The 
difference  between  these  two  signals  generates  an  error  voltage  that  drives 
a  voltage -controlled  oscillator  (VCO)  with  which  the  feedback  shift  register  is 

synchronized,  the reby  t racking  the  incoming  random  sequence  (Fig.  A-6). 

The  on-time  pseudo  random  sequence  is  mixed  with  the  incoming 
signal  to  reconstruct  the  carrier  signal.  The  phase  lock  loop  (which  is 
actually  a  bistable  Costas  loop)  tracks  this  satellite -transmitted  carrier. 
Binary  data  is  added  modulo -2  to  the  P  signal  pseudo  random  sequence  at 
a  rate  of  50  bps.  Since  only  the  pseudo  random  sequence  is  removed  from 
the  incoming  signal,  the  50-bps  data  sequence  still  remains  on  the  signal, 
but  the  single -sided  bandwidth  of  this  signal  is  now  only  about  50  Hz.  In 
addition  to  maintaining  phase  lock  on  the  carrier  signal,  the  Costas  loop 
(Fig.  A-7)  also  strips  off  the  data  remaining  on  the  signal. 
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Fig.  A -4.  Simplified  Diagram  of  Generic  Receiver 


Fig.  A -5.  Generic  Pseudo’  Random  Noise  Receiver 
Functional  Block  Diagram 
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Tracking  the  incoming  pseudo  random  sequence  allows  the  receiver 
to  make  a  range  measurement  to  the  satellite.  Tracking  the  carrier  signal, 
however,  does  not  provide  absolute  ranging  measurements;  only  changes 
in  the  range  can  be  measured  from  carrier  tracking.  Since  the  chipping 
rate  of  the  pseudo  random  sequence  corresponds  to  about  30  m  and  since 
the  wavelength  of  the  carrier  frequency  is  about  19  cm,  there  is  a  two-order - 
of-magnitude  difference  in  the  quantization  of  the  code  and  carrier  tracking 
loops,  and  the  accuracy  of  the  phase  lock  loop  is  much  greater  than  that  of 
the  code  loop. 
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APPENDIX  B 


NAVIGATION  EQUATIONS 


Figure  B-1  illustrates  an  earth -centered  inertial  coordinate  system. 
At  zero  time  the  x-axis  passes  through  the  intersection  of  the  equator  and 
prime  meridian,  the  z-axis  passes  through  the  North  Pole,  and  the  y-axis 
completes  the  right-handed  coordinate  system.  Because  of  earth  rotation, 
the  X  and  y  coordinates  change  in  longitude  about  15  deg  per  hour.  Shown  on 
the  figure  are  the  user  position  (x,  y,  and  z)  and  the  position  of  satellite  1 

(Xj,  y^,  and  z^).  The  range  distance  between  the  user  and  satellite  1  is 
shown  as  . 

The  basic  equations  using  four  satellites  are 

J(X  -  Xj)^  +  (y  -  y^)^  +  (z  .  +  T  = 

J(x  -  x^)^  +  (y  -  y^)^  +  (z  -  z^)^  +  =  ^2 

J{x  -  x^)^  +  (y  -  y^)^  +  (z-  +  T  =  R^ 

J(x  -  x^)^  +  (y  -  y^)^  +  (z  -  Z^)^  +  T  =  R^ 

where  x,  y,  z,  and  T  are  user  position  and  clock  bias  (unknowns),  x.,  y., 
and  z.  are  the  ith  satellite  position,  i  =  1,  4  (known),  and  R.  is  the  pseudo 
range  measurement  to  the  ith  satellite.  Here  the  quantities,  R^,  R^,  R^, 
and  R^  are  'pseudo  ’  ranges  in  that  they  are  the  sum  of  the  actual  range 
displacements  plus  the  offset  due  to  user  time  error.  For  convenience, 
units  have  been  selected  such  that  the  velocity  of  light  is  unity.  Roughly 
speaking,  if  displacement  is  measured  in  feet,  then  time  is  measured  in 
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Fig.  B-1.  Earth -Center ed  Coordinate 


nanoseconds  since  the  velocity  of  light  is  approximately  10^  ft/sec.  In  the 
equations  shown  here  the  four  pseudo  ranges  are  the  measured  quantities. 
The  satellite  positions  are  known,  and  the  four  unknowns  are  user  position 
and  the  user  clock  error.  It  should  be  emphasized  that  while  precision 
atomic  frequency  standards  are  used  in  the  satellites  and  the  monitor  sta¬ 
tions,  there  is  no  requirement  for  Navstar/GPS  users  to  have  a  precision 
clock.  Ordinary  quartz  crystal  frequency  standards  are  adequate  for  the 
user  since  he  is  continuously  computing  time  for  the  four  pseudo  range 
measurements . 

The  above  equations  are  nonlinear.  While  it  is  possible  to  solve 
these  equations  directly  as  they  are  shown,  user  equipments  without  ex¬ 
ception  employ  a  much  simpler  linearized  version  of  these  equations.  The 
basic  navigation  equations  can  be  linearized  by  employing  incremental 
relationships  as  follows. 

Let 

z^,  T^  be  nominal  (a  priori  best  estimate)  values 
of  X,  y,  z,  T 

Ay,  Az,  AT  be  the  corrections  to  these  nominal  values 
be  the  nominal  pseudo  range  measurements  from  the  ith  satellite 
ARj  be  the  difference  between  the  actual  and  nominal  measurements 
Therefore 

X 

y 

z 

T 

R. 

1 

R  . 
ni 
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=  X  +  Ax 
n 

=  y  +  Ay 

=  z  +  Az 
n 

=  T  +  AT 
n 

=  R  .  +  AR. 
ni  1 

=  -  y,)"  +  t  T_, 
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Substituting  the  incremental  expressions  into  the  basic  equations  yields 


+  Ax  -  x.)^  +  (y^  +  Ay  -  y. +  (z^  +  A^  -  z.)^ 

=  R  .  +  AR.  -  T  -  AT,  i  =  1,  4 
ni  in 


By  ignoring  second-order  error  terms,  these  equations  can  be  written  as 


J(x^  -  x.)^  +  (y^  -  y.)^  +  (z^  -  z/ 

(x  -  x.)Ax  +  (y  -  y*)^y  +  (z  -  z.)Az 
n  1  n _ i _ n _ i 


=  R  .  +  AR.  -  T  -  AT 
ni  in 


Substituting 


(x  -  X.)  (y^  -  y.) 

"  Ax  +  ^  Ay  + 


R  .  -  T 
ni  n 


R 


ni 


n 


(z  -  z.) 

R^^-^  Az  .  AT  =  AR. 


ni 


n 


The  above  four  equations  (i  =  1 ,  2,  3,  4)  are  the  linearized  equations 
that  relate  pseudo  range  measurements  to  the  desired  user  navigation  infor¬ 
mation  as  well  as  the  user's  clock  bias.  The  known  quantities  of  the  right- 
hand  side  of  the  equation  are  actually  incremental  pseudo  range  measure¬ 
ments.  They  are  the  differences  between  the  actual  measured  pseudo  ranges 
and  the  measurements  that  had  been  predicted  by  the  user's  computer  based 
on  the  knowledge  of  satellite  position  and  the  user's  most  current  estimate  of 
his  position  and  clock  bias.  The  quantities  to  be  computed.  Ax,  Ay,  Az,  and 
at,  are  corrections  that  the  user  will  make  to  his  current  estimate  of  posi¬ 
tion  and  clock  biases.  The  coefficients  of  these  quantities  on  the  left-hand 
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side  are  the  direction  cosines  of  the  line  of  sight  (LOS)  from  the  user  to 
the  satellite  as  projected  along  the  x,  y,  and  z  coordinates.  For  all  four 
equations  the  coefficient  in  front  of  AT  is  unity.  These  linearized  equations 
can  be  conveniently  expressed  in  matrix  notation  and  appear  as 


where  o  . ^  is  the  direction  cosine  of  the  angle  between  the  range  to  the  ith 
satellite  and  the  jth  coordinate. 


By  the  use  of  matrix  notation,  the  above  equations  can  be  expressed 
very  compactly  as  follows.  Let 

r  =  the  four -element  pseudo  range  measurement  difference  vector 
X  =  the  user  position  and  time  correction  vector 
A  =  the  4x4  solution  matrix 
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X  =  [Ax  Ay  Az  AT]^ 
r  =  [ARj  AR^  AR^  AR^]^ 

Therefore 

Ax  =  r  or  x  =  A”  ^r 

The  last  equation  presented  compactly  expresses  the  relationship 
between  pseudo  range  measurements  and  user  position  and  clock  bias.  Since 
this  relationship  is  linear,  it  can  be  used  to  express  the  relationship  between 
the  errors  in  pseudo  range  measurement  and  the  user  quantities.  This  re¬ 
lationship  is  therefore 


where  c  represents  the  pseudo  range  measurement  errors  and  e  ,  the  cor- 
.  ^ 
responding  errors  in  user  position  and  clock  bias. 

Let  us  now  consider  the  covariance  matrix  of  the  expected  errors  in 
pseudo  range  measurements  and  the  covariance  matrix  of  the  user  quantities. 
The  first  covariance  measurement  is  a  4  X  4  array  composed  of  the  expected 
values  of  the  squares  and  products  of  the  errors  in  the  pseudo  range  mea¬ 
surements.  The  diagonal  terms  in  the  matrix,  namely  the  squares  of  the 
expected  errors,  are  the  variances,  i.e.,  the  squares  of  the  expected  1  o 
values  of  the  pseudo  range  measurement  errors.  The  off-diagonal  terms  are 
the  covariance  between  the  pseudo  range  measurements  and  reflect  the  cor¬ 
relations  to  be  expected  in  these  measurements.  Likewise  the  covariance 
matrix  for  the  user  quantities  is  composed  of  the  expected  values  of  the 
squares  and  products  of  the  errors  in  the  user  quantities.  The  diagonal 
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terms  are  the  variance  or  the  squares  of  the  1  errors  in  user  position  and 
time,  while  the  off-diagonal  terms  reflect  the  correlations  in  these  errors. 
These  covariance  matrices  are  given  by 


COV(r)  = 

COV(x)  = 


where  the  symbol  e|  |  designates  "expected  value"  of  the  quantity  inside  the 
braces. 

Upon  substitution,  the  matrix  relationship  between  the  two  covariance 
matrices  becomes 

COV(x)  =  A"^COV(r)A“^ 

An  alternate  formulation  for  this  relationship  based  on  a  straight¬ 
forward  matrix  algebra  manipulation  is 

COV(x)  =  [A^COV(r)"^Aj"^ 

Note  that  the  relationship  between  the  pseudo  range  measurement  errors  and 
the  user's  position  and  clock  bias  errors  is  a  function  only  of  the  solution 
matrix  A,  which  in  turn  is  a  function  only  of  the  direction  cosines  of  the 
LOSs  from  the  user  to  the  satellites  along  whichever  coordinate  system  is 
used.  In  other  words,  the  error  relationships  are  functions  only  of  satellite 
geometry.  An  important  consideration  in  the  proper  use  of  Navstar/GPS  is 
the  fact  that  the  four  satellites  being  used  must  possess  "good  geometric 
properties.  (A  "good"  situation  is  one  in  which,  because  of  satellite 
geometry,  a  given  level  of  error  in  the  pseudo  range  measurements  results 
in  small  user  errors.  )  This  leads  to  the  concept  of  geometric  dilution  of 
precision  <GDOP),  a  measure  of  how  satellite  geometry  degrades  accuracy. 
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The  following  assumption  regarding  pseudo  range  measurement  errors 
provides  a  method  of  quantitatively  determining  whether  a  particular  four- 
satellite  geometry  is  good  or  bad.  Let  each  individual  pseudo  range  mea¬ 
surement  have  an  error  (la)  of  unity,  where  the  expected  mean  is  zero  and 
the  correlation  of  errors  between  satellites  is  also  zero.  With  these  as¬ 
sumptions  the  covariance  matrix  for  the  errors  in  the  pseudo  range  mea¬ 
surements  becomes  a  4  x  4  unity  matrix.  Thus,  for  this  case,  the  cova¬ 
riance  matrix  for  user  position  and  clock  bias  errors  is  given  by 

COV(x)  =  (A^A)'^ 


GDOP  is  defined  as  the  square  root  of  the  trace  of  COV(x)  when  COV(r) 
is  an  identity  matrix;  therefore 


Some  properties  of  this  quantity  can  be  sunamarized  as  follows: 

a.  GDOP  is,  in  effect,  the  amplification  factor  of  pseudo 
range  measurement  errors  into  user  errors  due  to  the 
effect  of  satellite  geometry. 

b.  GDOP  is  independent  of  the  coordinate  system  employed. 

c.  GDOP  is  a  criterion  for  designing  satellite  constellations. 

d.  GDOP  is  a  means  for  user  selection  of  the  four  best 
satellites  from  those  thar  are  visible. 

By  letting  V  ,  V  ,  V  ,  V.^  be  the  variances  of  user  position  and  time,  we 
X  y  2  A 

have 


GDOP  =  v/v  +  V  +  V  + 

V  X  y  z  T 
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As  an  alternative  to  GDOP  as  a  criterion  for  selecting  satellites  or 
evaluating  satellite  constellations,  only  some  of  the  variances  of  user  posi¬ 
tion  and  time  might  be  used.  These  are  defined  as  follows; 

PDOP  The  square  root  of  the  sum  of  the  squares  of  the 

three  components  of  position  error 

HDOP  The  square  root  of  the  sum  of  the  squares  of  the 

horizontal  components  of  position  error 

VDOP  The  altitude  error 

Note:  PDOP^  =  HDOP^  +  VDOP^ 

TDOP  The  error  in  the  user  clock  bias  multiplied  by  the 

velocity  of  light 

Note:  GDOP^  =  PDOP^  +  TDOP^ 

The  alternative  criterion  most  frequently  used  is  the  position  dilution 
of  precision  (PDOP).  PDOP  is  also  invariant  with  the  coordinate  system  and 
is  used  because  the  most  important  consideration  in  any  navigation  system  is 
position  accuracy;  knowing  time  is  generally  a  secondary  by-product.  Another 
alternative  is  the  horizontal  dilution  of  precision  (HDOP),  which  is  most 
meaningful  for  users  who  are  using  the  system  primarily  to  obtain  horizontal 
position. 
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APPENDIX  C 


BASELINE  ORBIT  CONFIGURATION 

The  space  segment  configuration  for  the  Phase  III  fully  operational 
system  is  summarized  as: 

a.  24  satellites 

b.  12 -hr  circular  orbits 

c.  5 5 -deg  inclination 

d.  Eight  satellites  in  each  of  three  orbit  planes 

e.  Satellites  in  a  plane  uniformly  spaced  45  deg  apart 

f.  Orbit  planes  uniformly  spaced  120  deg  apart  in  longitude 

g.  Satellites  staggered  15  deg  from  one  plane  to  another 

Table  C-1  shows  the  statistics  on  the  number  of  satellites  that  are 
visible  under  the  assumption  that  they  can  be  seen  above  either  5  or  10  deg 
in  elevation.  These  results  are  the  average  for  all  locations  on  earth  and 
for  all  times  during  the  day.  In  general,  eight  or  nine  satellites  can  be  seen 
at  a  given  time  above  6  deg,  and  six,  seven,  or  eight  can  be  seen  above 
10  deg.  Since  only  four  satellites  are  needed  for  navigation,  users  generally 
have  twice  as  many  available  from  which  to  select  the  best  four. 

The  geometric  performance  of  the  baseline  orbit  configuration  is 
shown  in  Fig.  C-1.  This  figure  shows  the  cumulative  probability  of  the 
position  dilution  of  precision  (PDOP)  for  three-dimensional  navigation 
using  four  satellites.  A  55-deg  inclination  angle  has  been  selected  to  op¬ 
timize  performance.  The  orbit  planes  intersect  each  other  at  about  90  deg 
and  form  eight  equal  octants. 


Table  C-1.  Navstar  Phase  III:  3x8  Distribution  of 
Global  Satellite  Visibility 


No .  of 
Satellites 
Visible 

5-deg  Elevation  Mask  Angle 

10-deg  Elevation  Mask  Angle 

Probability, 

percent 

Cunaulative 

Probability 

Probability, 

percent 

Cumulative 

Probability 

10 

4.  5 

4.  5 

0 

0 

9 

38.  5 

43.  0 

8.0 

8.  0 

8 

34.4 

77.4 

37.0 

45.  0 

7 

9.6 

87.  0 

23.6 

68.  6 

6 

13.  0 

100 

27.4 

96.  0 

5 

0 

100 

3.8 

99.8 

4 

0 

100 

0.2 

100 

42 


100 


I 
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Fig.  C-1.  Cumulative  Probability  of  PDOP 
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APPENDIX  D 

DOPPLER  NAVIGATION  WITH  NAVSTAR 

Navstar/GPS  performs  highly  accurate,  three-dimensional,  instanta¬ 
neous,  and  continuous  navigation  by  obtaining  pseudo  range  measurements 
from  four  satellites.  Theoretically,  pseudo  range  rate  measurements  could 
be  used  instead  of  pseudo  range,  and  this  idea  (using  integrated  pseudo  range 
rate,  i.e.,  measuring  the  Doppler  effect  on  the  transmission  from  the 
satellite  to  the  user  via  the  fully  operational  Phase  III  GPS)  has  been  investi¬ 
gated.  The  user  simiiltaneously  makes  an  integrated  Doppler  measurement 
from  four  satellites  over  an  interval  of  time  T.  (For  the  purpose  of  this 
analysis  it  is  assumed  that  the  time  intervals  of  the  measurements  from  the 
satellites  are  all  the  same  and  that  they  are  performed  simultaneously.  ) 

The  user  has  then  acquired  four  measurements  of  range  change  over  the 
interval  T,  and  from  these  measurements  he  can  compute  the  three  com¬ 
ponents  of  position  and  the  frequency  offset  of  the  local  oscillator  from 
nominal.  Obtaining  this  frequency  offset  is  analogous  to  the  user's  deter¬ 
mining  his  clock  bias  in  the  normal  case  of  making  pseudo  range  measure¬ 
ments  . 

For  the  above  approach  to  Doppler  navigation  the  position  dilution  of 
precision  (PDOP)  can  be  defined  as  follows;  Consider  that  the  four  range 
change  measurements  are  made  with  an  error  (lo)  of  unity  and  that  these  four 
measurement  errors  are  uncorrelated.  For  these  errors  there  will  be  a 
corresponding  three-dimensional  user  navigation  error  (la)  in  computing 
the  position:  PDOP  is  defined  as  this  three-dimensional  position  error  and 
is  the  ratio  of  user  error  to  measurement  error.  This  ratio  is  a  function 
of  the  relative  geometry  of  the  four  satellites  selected  to  perform  Doppler 
navigation.  To  obtain  the  best  possible  accuracy,  the  user  would  select  the 
set  of  four  satellites,  out  of  all  those  visible,  that  provide  the  minimum  PDOP. 
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Using  the  above  criteria  for  selecting  the  best  four  satellites,  a  com¬ 
puter  program  has  been  used  to  determine  statistically  the  PDOP  values 
that  can  be  expected  when  Doppler  navigation  is  employed.  The  results 
show  that  for  all  measurement  intervals  (T)  up  to  about  one  hour,  the  PDOP 
varies  inversely  with  T.  This  being  the  case,  a  PDOP  coefficient  can  be 
defined  by  multiplying  PDOP  by  the  measurement  interval.  In  this  manner 
the  statistical  results  can  be  condensed  and  are  readily  presented  for  the 
case  of  all  measurement  intervals.  The  cumulative  probability  distribution 
of  Kj,  which  is  PDOP  multiplied  by  the  measurement  interval  in  seconds,  is 
shown  in  Fig,  D-1.  The  PDOP  coefficient  ( l<r)  is  1.27X10”^,  and  PDOP  (la) 
is  therefore  given  by 


PDOP  =  1-27  X  10 

T 

be  the  error  (la)  in  the  integrated  pseudo  range  rate  measure¬ 
ment.  The  three-dimensional  user  position  error  is  given  by 

1.27  X  10^  X  CT 
IT  -  _ _ m 


Note  that  divided  by  T  is  the  average  pseudo  range  rate  error. 

To  summarize  the  statistical  results,  user  position  error  is  directly 
proportional  to  the  pseudo  range  rate  measurement  error ;  the  constant  of 
proportionality  is  1.27  x  10  sec.  For  example,  a  completely  stationary 
user  could  obtain  an  accurate  position  fix  in  2  min  if  he  measured  the  in¬ 
tegrated  pseudo  range  rate  by  tracking  the  1.  575-GHz  carrier  frequency. 
Assuming  that  this  measurement  is  made  with  an  accuracy  of  11  cm  over 
120  sec,  this  corresponds  to  a  measurement  error  of  about  0.09  cm/sec, 
ind  the  corresponding  user  position  error  will  therefore  be  about  12  m. 
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